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Abstract 

The Shape Calculus is a bio-inspired calculus for describing 3D 
shapes moving in a space. A shape forms a 3D process when combined 
with a behaviour. Behaviours are specified with a timed CCS-like pro- 
cess algebra using a notion of channel that models naturally binding 
sites on the surface of shapes. Processes can represent molecules or 
other mobile objects and can be part of networks of processes that 
move simultaneously and interact in a given geometrical space. The 
calculus embeds collision detection and response, binding of compati- 
ble 3D processes and splitting of previously established bonds. In this 
work the full formal timed operational semantics of the calculus is pro- 
vided, together with examples that illustrate the use of the calculus in 
a well-known biological scenario. Moreover, a result of well-formedness 
about the evolution of a given network of well-formed 3D processes is 
proved. 

1 Introduction 

In a near future, systems biology will profoundly affect health care and 
medical science. One aim is to design and test "in-silico" drugs giving rise to 
individualized medicines that will take into account physiology and genetic 
profiles |14| . The advantages of performing in-silico experiments by simulat- 
ing a model, instead of arranging expensive in- vivo or in- vitro experiments, 
are evident. But of course the models should be as faithful as possible to 
the real system. 

Towards the improving of the faithfulness of languages and models pro- 
posed in literature in the field of systems biology, the Shape Calculus [SOU] 
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was proposed as a very rich language to describe mainly, but not only bi- 
ological phenomena. The main new characteristics of this calculus are that 
it is spatial - with a geometric notion of a 3D space - and it is shape-based, 
i.e. entities have geometric simple or complex shapes that are related to 
their "functions", that is to say, in the context of formal calculi, the possi- 
ble interactions that can occur among the entities (called 3D processes, in 
the same context). Note that there are some other approaches that use a 
geometric notion of space [21 [161 H9 ; but they do not fully exploit the poten- 
tiality of geometry: they only consider positions that are centers of percep- 
tion/communication spheres or discretized the space using grids. Moreover, 
differently to classical mathematical models for biological systems, the Shape 
Calculus is individual-based, that is to say it considers autonomous entities 
that interact with others in order to give rise to an emerging behaviour of the 
whole system. However, this characteristic is also present in other languages 
and models proposed in literature |18| \19\ [TT1 [71 [12] . Parallel to the intro- 
duction of the Shape Calculus, a simulation environment called BloSHAPE 
|20[ [9] has been proposed. This environment embodies all the concepts and 
features of the Shape Calculus, plus others, more detailed, characteristics 
that are related to the implementation of simulations in different biological 
case studies EJ HDJ- The Shape Calculus can be considered the formal 
core of BloSHAPE. On the Shape Calculus, formal verification techniques 
can be applied, while in-silico experiments (simulations) can be performed 
in BloSHAPE. These two approaches can be used successfully in an inde- 
pendent way, but they can also be combined to address complex biological 
case studies. 

In [5] , the language of the Shape Calculus was introduced with the main 
aim of gently and incrementally present all its features and their relative 
semantics. The motivations behind the type and nature of the calculus op- 
erators were discussed and a great variety of scenarios in which the calculus 
may be used effectively were described. In [3] the timed operational seman- 
tics of the calculus was introduced and a well-formedness property of the 
dynamics of the calculus was proved. 

This paper is an extended version of [4j in which the original concepts 
of the calculus are fully recalled and in which further examples are pro- 
vided in order to explain better the technical details of its timed operational 
semantics. 

We present the full timed operational semantics of the Shape Calculus 
and we explain with proper running examples all the technical points that 
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Figure 1 : An example of binding and subsequent weak-split of two 3D pro- 
cesses. 

need particular attention and further explanation. Moreover, we define a 
concept of well-formedness, starting from shapes and porting it to more 
complex calculus objects such as 3D processes and, ultimately, networks 
of 3D processes. In the calculus, only well-formed objects are considered. 
Well-formedness is a standard concept used to avoid strange and unwanted 
situations in which a term can be legally written by syntax rules, but that 
semantically corresponds to a contradictory situation, for instance, in our 
composed shape whose pieces move in different directions. We prove 
that a given well-formed network of 3D processes always evolves into a well- 
formed network of 3D processes, that is to say, no temporal and spatial 
inconsistencies are introduced by the dynamics of the calculus. 

The paper is organized as follows. Section [2] recalls the main concepts 
of the calculus and focuses on weak and strong split operations whose se- 
mantics are slightly more complicated to define than that of other operators. 
Section [3] introduces formally 3D shapes, shape composition, movement, col- 
lision detection and collision response. Section [4] defines behaviours and 3D 
processes giving them full semantics. Section [5] puts all the pieces together 
and specifies precisely networks of 3D processes and a general result of dy- 
namic well-formedness is proven. Finally, Section [6] concludes with ongoing 
and future work directions. For the sake of readability all the proofs have 
been moved in Appendices A, B and C. 
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2 Recall of Main Concepts 



In this section we recall the main concepts the calculus introduced in [5] and 
we focus in particular to the weak and strong split operations, in order to 
make them clear from the beginning and, thus, then smoothly define their 
semantics, which requires some particular technical expedients. 

The general idea of the Shape Calculus is to consider a 3D space in 
which several shapes reside, move and interact. While time flows, shapes 
move according to their velocities, that can change over time both due to a 
general motion law - for instance as in a gravitational or in an electromag- 
netic field, or as in Brownian motion - and due to collisions that can occur 
between two or more shapes. Collisions can result in a bounce, that is to 
say elastic collisions. Instead, as it often happens in biological scenarios, 
colliding objects can bind and become a new compound object moving in 
a different way and possibly having a different behaviour. In this case we 
speak of inelastic collisions since they are treated with the physical law for 
that kind of collision. 

The Shape Calculus can be used to represent a lot of scenarios at differ- 
ent scales in different fields [5j. However, in this work we use a well-known 
biological scenario in order to introduce simple running examples that il- 
lustrate the semantics of the calculus operators. We consider biochemical 
reactions occurring inside a cell. Every species of involved molecules has a 
specific shape and we know from biology that the functions of a molecule 
are tightly related to its shape. For instance, in enzymatic reactions the 
functional sites that are active in the enzyme structure, at a given time, 
determine which substrate (one or two metabolites) can bind the enzyme 
and proceed to the catalyzed reaction. 

Fig. [T] shows a (2D, for simplicity) possible representation of an enzy- 
matic reaction. The larger object represents an enzyme with a shape So 
and a behaviour Bq; So and Bo together constitute a 3D process So[Bq]. 
The 3D process Si[-Bi] represents a metabolite that is close to the given 
enzyme. Note that portions of the shape surfaces are highlighted: they are 
the channels (an extension of the notion of channels of CCS [17J) that the 
corresponding 3D processes exhibit to the environment. Each channel has 
a name and an active surface. For instance, (a, X) is a channel of type a 
on the active site X of So [Bo] ■ The enzyme in Fig. [I] has two channels and 
its behaviour can be specified as: Bq = {b,Y).B' Q + {a,X).B r Q. The plus 
operator represents a non-deterministic choice between two possible com- 
munications on the channels. This non-determinism is resolved during the 
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Figure 2: An example of complex formation and subsequent strong-split. 



evolution of the system depending on which 3D processes will collide with 
the enzyme and where. 

Following the evolution proposed in the figure, after some time t elapsed 
(represented by the transition A) and after a detection and resolution of an 
inelastic collision (transition —>•$), we get one compound process represented 
by the term So[B' ](b, W)Si[B' 1 \. The bond is established on the surface of 
contact W and the name b records the type of channels that bound. Note 
that communication, i.e. binding, can only happen if there is a collision be- 
tween two 3D processes that expose compatible channels (name and co-name 
a la CCS) on their surface of contact. If the channels were not compatible, 
the collision would have been treated as elastic and the two 3D processes 
would have simply bounced. By letting B[ = (c,Yi).B" we allow the com- 
ponent of shape S\ to open a new channel and to bind with other colling 
processes. 

The third stage of Fig. [T] represents a possible evolution of an enzyme 
binding with a substrate; it can happen that, for some reason, the bond 
is loose and the two molecules are free again. To model this kind of be- 
haviour we introduce a special kind of not urgent split operation, called 
weak-split operation, that can be delayed of an unspecified time. This is an- 
other source of non-determinism in the calculus. Fig. [2] shows what happens 
when another substrate - that we represent as the 3D process S2 [B2] with a 
channel {a, X\) - binds with the compound process So[B' ]{b, W)Si[B[] on 
the common surface W\. In the terminology of biochemical reactions, a final 
complex has been formed. As a consequence, the reaction must proceed and 
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the products must be released. In our calculus, reactions are represented 
by strong-split operations. Differently from weak-split operations, this kind 
of split operations cannot be delayed and must occur as soon as they are 
enabled, i.e. when all the involved components can release all the bonds. In 
this example the involved components are Sq[B' ], Si[B[] and S^f-B^]; the 
set of bonds to be split is L = {(b, W), (a, Wi)}. 

3 3D Shapes 

We start by introducing three dimensional shapes as terms of a suitable 
language, allowing simpler shapes to bind and form more complex shapes. 
From now on we consider assigned a global coordinate system in a three 
dimensional space represented by M 3 . Let P,V = M 3 be the sets of posi- 
tions and velocities, respectively, in this coordinate system. Throughout the 
paper, we assume relative coordinate systems that will always be w.r.t. a 
certain shape S, i.e. the origin of the relative system is the reference point 
p of S. We refer to this relative system as the local coordinate system of 
shape S. If p £ P is a position expressed in global coordinates, and V C P 
is a set of points expressed w.r.t. a local coordinate system whose origin is 
p, we define global(V, p) = V + p = {x + p | x G V} to be the set V w.r.t. 
the global coordinates. Using local coordinate systems we can express parts 
of a given shape, such as faces and vertexes, independently from its actual 
global position. 

Definition 1 (Basic Shapes) A basic shape a is a tuple {V, m, p, v) where 
V C P is a convex polyhedron (e.g. a sphere, a cone, a cylinder, ete.J^J that 
represents the set of shape points; m € M + , p G P and v G V are, resp., the 
mass, the centre of mass^and the velocity of a. All possible basic shapes are 
ranged over by a, a', ■■■ . We also define the boundary B(cr) of a to be the 
subset of points of V that are on the surface of 

Note that basic shapes are very simple convex shapes. They can be repre- 
sented by means of suitable and efficient data structures and are handled by 
the most popular algorithms for motion simulation, collision detection and 
response [13J. Three dimensional shapes of any form can be approximated 

2 From a syntactical representation point of view, we assume that V is finitely repre- 
sented by a suitable data structure, such as a formula or a set of vertices. 

3 We actually need only a reference point. Thus, any other point in V can be chosen. 
4 Note that we consider only closed shapes, i.e. they contain their boundary. 
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- with arbitrary precision - by composing basic shapes: the composition of 
two shapes corresponds to the construction of a third shape by "gluing the 
two components on a common surface. This concept is generalized by the 
following definition. 

Definition 2 (3D shapes) The set S of the 3D shapes, ranged over by 
S, S', ■ ■ ■ , is generated by the grammar: 

S ::=a\ S (X) S 

where a is a basic shape and X C P. If S = a = (V, m, p,v), we define 
V{S) = V , m(S) = m, 7Z(S) = p, v(S) = {v} to be, resp., the set of points, 
the mass, the reference point and the velocity of S. If S = S\ {X) S 2 is a 
compound shape, then: V(S) = V(Si) U V{S 2 ), m(S) = m(Si) + m{S2), 
K{S) = (m(5i) • n(S!) + m(S 2 ) ■ W(5 a ))/(m(fi'i) + m{S 2 )f\ and v(5) = 
v(S*i) U v(S < 2). Finally, the boundary of S is defined to be the set B(S) = 
(B(Si) U B(S 2 )) \ {x £ P I x k interior ofV{S 1 (X) S 2 )}, where a point 
x € V C P is said to be interior if there exists an open ball with centre x 
which is completely contained in V. 

In this paper we only consider shapes that are well-formed according to 
the following definition. 

Definition 3 (Well-formed shapes) Each basic shape a is well-formed. 
A compound shape Si (X) S 2 is well-formed iff: 

1. both S\ and S 2 are well-formed; 

2. the set X = V(S\) n V{S 2 ) is non-empty and equal to B{S\) n B(S 2 ); 

3. v(S% (X) S 2 ) is a singleton {v]|^] where v = v(5i) = v(S*2). 

Below we say that two shapes S\ and S 2 interpenetrate each other if there 
exists a point x that is interior of both V(S\) and V^S^. In other terms, 
they interpenetrate iff the set X = V(S\) Pi V{S 2 ) 7^ is not a subset of 
B(S%) n B(S 2 ). If S = Si (X) S 2 is well-formed, X is said to be the surface 
of contact between S± and S 2 ; moreover, each x G X is a point of contact. 

Again for simplicity, we use the centre of mass as the reference point. Any other point 
can also be chosen. 

6 With abuse of notation, throughout the paper, we write v(5") also to refer to the 
element v of the singleton {v} as this is not ambiguous when S is well-formed. 
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Condition (|2]) guarantees that well-formed compound shapes touch but 
do not interpenetrate; the surface of contact X is always on the boundary of 
both S\ and S 2 . It can be a single point, a segment or a surface, depending 
on where the two shapes are touching without interpenetrating. Most of 
the time X is a (subset of) a feature of the basic shapes composing the 
3D shape, i.e., a face, an edge or a vertex. Condition (|3j imposes that all 
the shapes forming a compound shape have the same velocity; thus, the 
compound shape moves as a unique body. 

A compound 3D shape S can be represented in a number of different 
ways by rearranging its basic shapes and surfaces of contact. All these pos- 
sible representation are 'equivalent' w.r.t. the structural congruence defined 
below. 

Definition 4 (Structural Congruence of 3D Shapes) The structural 
congruence relation over S, denoted by =s, is the smallest relation that sat- 
isfies the following rules: 

1. S l (X)S 2 = s S 2 (X)S i; 

2. {Si (X) S 2 ) (Y) S 3 = s Si {X} (S 2 (Y) 5 3 ) provided that the surface of 
contact Y C B(S 2 ) n B(S$). 

The next example shows how the particular features of Shape Calculus 
can be used to construct a model in a well-known biological scenario. 

Example 1 (A Biological Example) The glycolysis pathway is part of 
the process by which individual cells produce and consume nutrient molecules. 
It consists of ten sequential reactions, all catalyzed by a specific enzyme. We 
focus, in this example, on the first reaction that can be described as 

glucose, ATP glucose-6-phosphate, ADP, H + 

where an ATP is consumed and a molecule of glucose (GLC) is phosphory- 
lated to glucose 6-phosphate ( G6P), releasing an ADP molecule and a positive 
hydrogen ion (Hydron). 

The enzyme catalyzing this first reaction is Hexokinase (HEX). GLC, 
G6P, ATP, ADP and H + are metabolites. Both enzymes and metabolites 
are autonomous cellular entities that continuously move within the cyto- 
plasm. The transformation of a metabolite into the one that follows in the 
"pipeline" of the pathway (in this case, GLC into G6P) depends on the meet- 
ing (collision in binding sites) of the right enzyme (in this example HEX) 
with the right metabolites, in this example GLC and ATP. The order of these 
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Figure 3: A network of 3D processes for describing the first reaction of the 
glycolysis pathway. 



bindings does not matter. After this binding the reaction takes place and the 
product^ are released in the cytoplasm. A special case occurs when the en- 
zyme has bound one metabolite and an environmental event makes it release 
the metabolite and not proceed to the completion of the reaction. 

We model the shape of HEX, which we denote with S^, by a polyhedron 
approximating its real shape and mass, available at public databases (e.g. 
Fig. [3] shows a network of 3D processes in which there are two hexokinases 
and some GLC, G6P, ATP and ADP 3D processes. Note that we use a 
unique kind of shape for GLC ad G6P, denoted by S g , and another unique 
kind of shape for ATP and ADP, denoted S a - They will be distinguished by 
their behaviours. 

Note that modelling of biochemical reactions with our calculus is very 
different from the usual ODE-based approach. This is because the Shape 
Calculus embeds concepts and features of a particle-based, individual-based 
and spatial geometric approach and we wanted to show, in a relatively simple 
and well-known scenario, how all these concepts and features can be used. 



3.1 Trajectories of Shapes 

The general idea of the Shape Calculus is to consider a three-dimensional 
space in which several shapes reside, move and interact. One of the choices 

7 In this example we neglect the hydron. 
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to be made is how the velocity of each shape changes over time. We believe 
that a continuous updating of the velocity - that would be a candidate 
for an 'as precise as possible approach of modelling - is not a convenient 
choice. The main reason is the well-known compromise between the benefits 
of approximation and the complexity of precision. Our choice, quite common 
also in computer graphics |13| . is to approximate a continuous trajectory of 
a shape with a polygonal chain, i.e. a piecewise linear curve in which each 
segment is the result of a movement with a constant velocity. The vertices 
of the chain corresponds to velocity updates. 

To this purpose we define a global parameter A G M + , called movement 
time step, that represents the maximum period of time after which the ve- 
locities of all shapes are updated. The quantification of A depends on the 
desired degree of approximation and also on other parameters connected to 



collision detection (see Section 3.2). In some situations, the time of up- 
dating can be shorter than A because, before that time, collisions between 
moving shapes can occur. These collisions must be resolved and the whole 
system must re-adapt itself to the new situation. The time domain T = Mq" 
is then divided into an infinite sequence of time steps U such that to = 
and ti < ti-i + A for all i > 0. An example in |5] (cf. Section 2) shows in 
more details how the timeline can be broken up into such time instants. 

In the calculus, the updating of velocities is performed by exploiting a 
function steer : T — > (S ^->- V) that describes how the velocity of all existing 
shapes (i.e. all shapes that are currently moving in the space) at each time 
t is changed. We assume that, at any given time instant t E T, steer t S is 
undefined iff shape S does not exist and, hence, its velocity has not to be 
changed. 

This approach provides us with the maximal flexibility for defining mo- 
tion. Static shapes are those shapes whose velocity is always zercj^j A 
gravity field can be simulated by updating the velocities according to the 
gravity acceleration. A Brownian motion can be simulated by choosing a 
random 3D direction and then defining the length of the vector w.r.t. the 
mass and/or the volume of the shape. In this paper, we do not consider 
movements due to rotations. However, this kind of movements can be easily 
added to our shapes by assigning an angular velocity and a moment of in- 
ertia to a shape and then by performing a compound motion of translation 
and rotation along the movement time step. 



8 To represent walls, we also need to assign an infinite value to the mass of these objects, 
otherwise they can be moved anyway due to collisions. 
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Let us define now some useful notation and properties. 

Definition 5 (Evolution of shapes over time) Let S G S and t G T; 

S + t, i.e. the shape S after t time units, is defined by induction on S: 

Basic: (V, m, p, v) + 1 = (V + (i • v), m, p + (t ■ v), v) 
Comp: (Si (X) S 2 ) + t= (Si +t)(X + t- v(S)) (S 2 + t) 

Definition 6 (Updating shape velocity) Let S G § anc? w € V. We 
define the shape S|[w|, i.e. S whose velocity is updated with w, as follows: 

Basic: (V, m, p, v) |[w| = (V, m, p, w) 

Comp: (Si (X) S 2 )|[w] = (Si|[w|) (X) (S 2 |[w]|) 

The following result comes directly from Def. [3] 

Proposition 1 Let S G S, t G T anc? w G V. If S is well-formed then S + t 
and S\[w]\ are well-formed. 

Our intent is to represent a lot of shapes moving simultaneously in space 
as described above. Inevitably, this scenario produces collisions between 
shapes when their trajectories encounter. 

3.2 Collision Detection 

There is a rich literature on collision detection systems, as this problem is 
fundamental in popular applications like computer games. Good introduc- 
tions to existing methods for efficient collision detection are available and 
we refer to Ericson |13] and references therein for a detailed treatment. 

For our purposes, it is sufficient to define an interface between our cal- 
culus and a typical collision detection system. We can then choose one of 
them according to their different characteristics, e.g. their applicability in 
large-scale environments or their robustness. It must be said, however, that 
the choice of the collision detection system may influence the kind of ba- 
sic (or compound) shapes we can use, as, for instance, some systems may 
require the use of only convex shapes to be more efficient^] 

Typically, a collision detection algorithm assumes to start in a situation 
in which shapes do not interpenetrate. Then it tries to move all the shapes 

9 The basic shapes that we consider in Definition [l] are typically accepted by most of 
the collision detection systems. 
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Figure 4: Some steps to determine the first time of contact between two 
shapes. 

of a littlef^l time step - that we have already introduced as the movement 
time step A - and check if interpenetrations occurred/^"] If so, it tries to 
consider only half of the original time step and repeat the interpentration 
check, i.e. it performs a binary search of the first time of contact t between 
two or more shapes, with some degree of approximation. Fig. [4] shows these 
steps. In case (a) the passage of the whole A results in an interpenetration. 
Then, in (b) the passage of A/2 is tried resulting into a non-contact. After 
some iterations the situation in (c) is reached. 

In addition to the first time of contact, a collision detection algorithm 
usually outputs the shapes that are colliding, i.e. are touching without in- 
terpenetrating after t, and some information about the surfaces or points of 
contact. We now define precisely what we expect to obtain from a collision 
detection system. 

Definition 7 (First time of contact) Let I be a non-empty finite set of 
indexes and let {Si}i^j be a set of well- formed shapes such that for alli,j G /, 
Si and Sj do not interpenetrate (see Def. [$p. The first time of contact of 
the shapes Si, denoted Ftoc({<Sj}j e i), is a number t £ T such that: 

10 The time step must be chosen little enough to avoid that any two shapes, at their 
possible maximum velocity, start in a non-interpenetration state, engage in an interpen- 
etration and exit from the interpenetration state in a single time step duration. 

11 Typically, the major efforts of optimization are concentrated in this step since the 
number of checks is, in the worst case, 0(N 2 ) - where N is the number of shapes in the 
space - but the shapes that are likely to collide are only those that are very close to each 
other. 
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1. for all t' E T, < t' < t and for all i, j £ I, Si + t' and Sj + t' do not 
interpenetrate; 

2. there exist i,j E /, with i ^ j, such that B(Si + t) DB(Sj +t) ^ 0, i.e., 
some shapes are touching at t; 

3. for all e > there exists 5, < 5 < e, and i,j G /, with i ^ j, such 
that Si + (t + 5) and Sj + (t + 5) interpenetrate, i.e., in t some shapes 
are touching and any further movement makes them to interpenetrate. 

Note that such a definition allows shapes that are touching without in- 
terpenetrating, and with velocities that do not make them to interpenetrate 
(e.g., the same velocity), to move without triggering a first time of contact. 
This will be useful in Section [5] when we split previously compound shapes. 
Indeed, after the split these shapes have the same velocity and, hence, do 
not affect the next first time of contact. 

Definition 8 (Collision information) Let {5j},; 6 / be a set of well-formed 
shapes and let t = Ftoc({5j}j g /) be their first time of contact. The set of 
colliding shapes after time t is denoted by colliding({S'j}j e /) CSxSx 
A tuple (Si,Sj,X) G colliding({5i} ie/ ) iff: 

1. V(Si + t)nV(Sj + t) is non-empty and is equal to B(Si + t)r\B(Sj+t); 

2. for all e > there exists 5, < 5 < e, such that Si + (t + 5) and 
Sj + (t + 5) interpenetrate. 

3.3 Collision Response 

In this section, we briefly discuss the problem of collisions response |15j . 
i.e. how collisions, once detected, can be resolved. We distinguish between 
elastic collisions (those in which there is no loss in kinetic energy) and per- 
fectly inelastic ones (in which kinetic energy is fully dissipated)^] After an 
elastic collision, the two shapes will proceed independently to each other 
but their velocities will be changed according to the laws for conservation 
of linear momentum and kinetic energy - Equations (l)-(2) in Def. [9] On 
the contrary, two shapes that collide inelastically will bind together and will 
move as a unique body whose velocity is given by the law for conservation 
of linear momentum only - Equation (3), in Def. [9] 

12 Other different kinds of collisions can be easily added to the calculus provided that 
the corresponding collision response laws are given. 
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Definition 9 (Collisions) Let S\,S2 £ § and let X C F be a surface of 
contact between them. If X is neither an edge nor a vertex of Si, the ve- 
locities wi and w 2 of these shapes after an elastic collision in X are given 
by: 

(1) wi=v(5i) 7TrT- n ( 2 ) w 2 = v(5 2 ) + , A > • n 

m[bi) m{b2) 

where n is the normal of contact away from X C B(S\), i.e. the unit vector 
perpendicular to the face of S\ that contains X, and 

x _ 2 m(Si)m(S 2 ) v(ffi) • n - v(6 , 2 ) • n 
m(Si) + m(5 2 ) nn 

If X is either an edge or a vertex of Si, n is the normal of contact away from 
the shape S" 2 and velocities wi and w 2 are obtained by means of symmetric 
equations. In both cases, we write Si < — > e S2 to denote the pair of velocities 
(wi,w 2 ). If Si and S2 collide inelastically in the surface of contact X, they 

will bind tog 
is given by: 



will bind together as a unique body whose velocity (denoted with Si < — >i S2) 



(3 v = — — — • v(5i) + t^tt-v5 2 

m(6i)+m(6 2 ) m(6i)+m(A 2 ) 



4 3D Processes 



In this section we introduce the timed process algebra whose terms de- 
scribe the internal behaviour of our 3D shapes. This is a variation of 
Timed CCS |21| . where basic actions provide information about binding 
capability and splits of shape bonds. We use the following notation. Let 
A = {a, 6, • • • } be a countably infinite set of channels names (names, for 
short) and A = {a \ a £ A} its complementation; by convention we assume 
a = a for each name a. Elements in A = A U A are ranged over by a, ft, ■ ■ ■ . 

Binding capabilities are represented by channels, i.e. pairs (a, X) where 
a G A is a name and X is a surface of contact. Intuitively, a surface of 
contact is a portion of space (usually, a subset of the boundary of a given 
3D shape) where the channel itself is active and where binding with other 
processes are possible. Names introduce a notion of compatibility between 
channels useful to distinguish between elastic and inelastic collisions. If (3 = 
a and X n Y 7^ we say that the channels (a, X) and (f3, Y) are compatible, 
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written (a, X) ~ Y). Otherwise, {a, X) and (/3,Y) are incompatible 
which we write as (a,X) 76 (f3,Y). 

We also introduce two different kinds of actions that represent splits of 
shape bonds. In particular, we distinguish between weak-split actions of the 
form Lo(a, X) and strong-split actions of the form p(a,X). With an abuse 
of notation, two weak-split actions u)(a,X) and oj(f3,Y) (as similarly for 
the strong-split actions p(a, X) and p(/3,Y)) are compatible if so are the 
channels (a, X) and {(3,Y). We will see that a synchronization between a 
pair of compatible weak-split actions results in a weak-split operation, while 
synchronizations between multiple pairs of compatible strong-split actions 
correspond to a strong-split operation. These operations behave differently 
w.r.t. to time passing since the latter operation cannot time pass further, 
while the former one can be arbitrarily delayed. 

Let C be the set of all channels, u(C) = {u:(a, X) \ (a, X) G C} and 
p(C) = {p(a, X) | (a, X) G C} be the sets of weak-split actions and strong- 
split actions, resp. Our processes perform basic and atomic actions that 
belong to the set Act = C U oj(C) U p(C) whose elements are ranged over 
by /i, //,■■■ . We finally assume a countably infinite collection K. of process 
name or process constants. 

Definition 10 (Shape behaviours) 

The set of shape behaviours, denoted by B, is generated by the following 
grammar: 

B :: = nil | {a,X).B \ Lu(a,X).B \ p{L).B \ e(t).B \ B + B\ K 

where (a,X) G C, L C C (non-empty) whose elements are pairwise incom- 
patible (i.e. for each pair (a,X),(/3,Y) G L it is (a,X) (f3,Y)), t £ T 
and K is a process name in K.. 

A brief description of our operators now follows, nil is the Nil-behaviour, 
it can not perform any action but can let time pass without limits. A 
trailing nil will often be omitted, so e.g. we write (a, X).uj(a, X) to abbreviate 
(a,X).oj(a,X).n\\. (a,X).B and uj(a,X).B are (action-) prefixing known 
from CCS; they evolve in B by performing the actions (a,X) and uj(a,X), 
resp. {a,X).B exhibits a binding capability along the channel (a, X), while 
oj(a, X).B models the behaviour of a shape that, before evolving in B, wants 
to split a single bond established via the channel (a, X). p(L).B is the 
strong-split operator; it can evolve in B only after the execution of all strong- 
split actions p(a,X) with (a,X) G L. The delay-prefixing operator e(t).B 
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p e Cuw(C) 

Nil, Pref, Str+ 

t t t 

nil -w nil p.B ~» p.B p(L).B ~» p(L).B 

B 1 B{ B 2 & B' 2 t' >t 

SUMj DELr 



B 1 + B 2 ^ B[ + B' 2 e(t').B -$> e{t' -t).B 

B ~*^> B' 
K ■$* B' 



B ~~> B' d f 

Def< if K = B 



Table 1: Temporal behaviour of B's terms 

(see |21] ) introduces time delays in 3D processes; t S T is the amount of 
time that has to elapse before the idling time is over (see rules DELAYj and 
DELAY a in Tables [T] and [2| . Finally, B\ + B 2 models a non-deterministic 
choice between B\ and B 2 and K is a process definition. 

In the remainder of this paper, we use processes in B to define the internal 
behaviour of 3D shapes. For this reason, we assume that sites in binding 
capabilities, as well as in weak- and strong-split actions, are expressed w.r.t. 
the local coordinate system whose origin is the reference point of the shape 
where they are embedded in. 

Definition 11 (Operational semantics of shape behaviours) The 

SOS-rules that define the temporal transition relations w C (B x B) for 
t 6 T, that describe how shape behaviours evolve by letting time t pass, are 
provided in Table i| We write B B' if (B,B') G-w and B if there is 
B' € B such that (B, B') G~~». Similar conventions will apply later on. Rules 
in Table^ define the action transition relations Ac (B x B) for fi £ Act. 
These transitions describe which basic actions a shape behaviour can perform. 

Most of the rules in Table [T] are those provided in |21] , Rules Pref^ 
and StRj state that processes like {a,X).B, u(a, X).B and p{L).B can be 
arbitrarily delayed. The only rules in Table [2] worth noting are those defin- 
ing the functional behaviour of the strong-split operator. By Rules StRi 
and STR2, if (ct,X) € L then p{L).B can do a p(a, X)-action and evolve 
either in B (if L = {(a,X)}) or in p(L\{(a, X)}).B (otherwise). Rule Str 3 
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H € C U w(C) B A B' B x A B' 

PREF a DEL a SUM a 

/i.B A B e(0).B A^ 5i + B 2 A B' 



BAB' def L = {(a,X>} 

DEF a ifK = B STRi- 



L = {(a, X)} U L' L' ^ B B' 

Str 2 STR; 



^.S p(L').B p(L).B p(L).B> 



Table 2: Functional behaviour of B-terms 

is needed to handle arbitrarily nested terms, e.g. p({(a, X}).p({(b,Y}).B. 
Other rules are as expected. 

Now we are ready to define 3D processes, i.e. simple or compound shapes 
whose behaviour is expressed by a process in B. 

Definition 12 (3D processes) The set 3DP o/3D processes is generated 
by the following grammar: 

P ::= S[B] | P(a,X)P 

where S € 8, B G B, a £ A and X is a non-empty subset of P. The shape 
of each P G 3DP is defined by induction on P as follows: 

Basic: shape(S'[B]) = S 

Comp: shape(P (a, X) Q) = shape(P) (X) shape(Q) 

We also define v(P) = v(shape(P)) and B(P) = B(shape(P)). Below we of- 
ten write P AAj Q and P AA e Q as shorthand for shape(P) AAj shape(Q) 
and shape(P) < — > e shape(Q), resp. Finally, P|[v| is the 3D process we ob- 
tain by updating P 's velocity as follows: 

Baszc: (5[B])|[v]| = (S[v])[B] 

Comp: (P (a, X) Q)|[v]| = (P|[v|) (a, X) (Q|[v]|) 

We finally write steer tP to denote P|[steer tshape(P)]. We say that a 
basic process S[B] is well-formed iff the shape S is well-formed and, for each 
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I C P that occurs in B, global(X, TZ(S)) C B(S). A compound process 
P (a, X) Q is well-formed iff P and Q are well-formed, v(P) = v(Q) and 
the site X expressed w.r.t a global coordinate system is a non-empty subset 
of B{P) (~)B(Q). Note that this also means that the set V(P) nV(Q) is non- 
empty and equal to B{P) D B{Q). Later on in this paper we only consider 
well-formed processes. 

We can state the following proposition as an easy consequence of shapes 
and 3D processes well-formedness. 

Proposition 2 For each P E 3DP well-formed, shape(P) is well-formed. 

Let us model the molecules involved in the reaction of Example [T] as 3D 
processes. 

Example 2 (3D Processes for HEX, GLC and ATP) An Hexokinase 
molecule is modeled as HEX = f S/JHEX] where: 
HEX=(atp,X te ).HA+(glc,X /ls ).HG ; 

HA = w(atp J X ha ).HEX + e(t h ).(glc,X /l| ,).p({(atp,X ho ),(glc,y Aj ,)}).HEX, 
HG = w(glc, X hg )MEX + e(ifc,).(atp, X ha ).p({(atp, X ha ), (glc, Y hg )}).HEX, 

and Xh a ,Yhg are the surfaces of contact shown in Fig. [3| ATP = S a [ftTP] 
models an ATP molecule where: 

ATP = (^,X ah ).(e(t a ).p({(^,X ah )}).ADP + oj(^,X ah )AJP) 

and the surface of contact X a h is the whole boundary B(S a ).The process 
modelling a molecule of glucose is similar: GLC = S g [GLC] where 

GLC = (gk, X gh ).(e(t g ).p({@fc, X gh )}).G6P + w(gk, X gh ).GLC 
We leave unspecified the behaviours G6P and ADP. 

HEX has two channels (atp,Xh a ) and (glc, Yh g ) to bind, resp., with an 
ATP and a GLC molecule. By performing an action (atp, X^a), HEX evolves 
in HA. HA can perform either a weak-split action u;(atp, Xha) to come back 
to HEX, or can wait at most th units of time, perform (glc, Y^g) and then 
evolve in /)({(atp, X/,, a ), (glc, Y/ l9 )}).HEX. Now, two strong-split actions are 
enabled after which we come back to HEX. Notice that, after an action 
(glc, Y^g), HEX becomes HG that behaves symmetrically. 

An ATP molecule performs a (atp, X a h) -action, waits t r units of time, 
and then can release the bond established on the channel (atp, X a f- L ) - and 
thus return free as ATP - or can participate in the reaction and become an 
ADP. As we will see in Section^ the result is the split of the complex in the 
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three original shapes whose behaviours are HEX, ADP and G6P, resp. We 
omit the description of the behaviour of GLC since it is similar to that of 
ATP. 

We are ready to define the timed operational semantics of 3D processes. 

Definition 13 (Transitional semantics of 3D processes) Rules in Ta- 
ble [5] define the transition relations ic (3DP x 3DP) for t € T, and 
— >C 3DP x 3DP for fi € Act. Two 3D processes P and Q are said to 

be compatible, written P ~ Q, if P — — -> and Q - ' ^> for some compatible 
channels (ct,X) and (a,Y); otherwise, P and Q are incompatible that we 

denote with P 76 Q. Below, we often write P -jk and P ^ as shorthand for 

a(a,X) ty(a,X) . . 

P f — : — > and P f — ■ — >, resp., for any {a,X). 

Essentially, rules in Table [3] say that a 3D process inherits its functional 
and temporal behaviour from the IB-terms defining its internal behaviour. 
But now sites of binding capabilities and split actions are expressed w.r.t. 
a global coordinate system (see rules BASIC C and BASIC S ). For simplicity, 
we have omitted a rule defining which weak-split action a basic process can 
perform. This can be obtained from rule BASIC^ by replacing each pQ- 
action with a corresponding wQ-action. It is worth noting that, due to rule 
COMP a 2, some of the (a, y)-actions performed by P (by Q) can be prevented 
in P (a, X) Q since, due to binding a part (or all) of Y has became interior 
because it is covered by a piece of Q (of P respectively) the surface of contact 
Y ^ B(P (a, X) Q) and, hence the corresponding channel is no more active. 
We have also omitted symmetric rules for COMP„i and COMP a 2 for the 
actions of Q. 

The following proposition (see Appendix B for the proof) shows that 3D 
processes well-formedness is closed w.r.t. transitions and A. 

Proposition 3 Let P G 3DP well-formed. Either P Q or P A Q implies 
Q £ 3DP well-formed. 

At this stage a key observation is that the operational rules in Table [3] 
do not allow synchronization between components of compound process that 
proceed independently to each other. Consider, as an example, 

P = S p [p({^X p }).B p ], 
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Basic* 

S[B] - (S + t)[B'] 

p ^p> Q^Q* X' = X + (t ■ v(P)) 

COMP t 

P (a, X)Q^P' (a, X') Q' 



B ^% B' Y = global(X, K(S)) 

BASIC C -—77 

S[B] -^-4 S[B'} 
B^%B' Y = g\ob*\(X,n(S)) 

BASIC S — — 

S[B] ^-4 S[B'} 



COMP s - 



p P(°>Y\ pi 

P(a,X)Q^%P'(a,X)Q 
p ^»' y > ; p i 



COMP™ — 

P{a,X)Q^-U P> (a,X)Q 

P^Xp' Y CB(P(a,X)Q) 

COMPc (aY) 

P(a,I)Q^4P' (a,X)Q 



Table 3: Functional and temporal behaviour of 3DP-terms 
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Q = S q [p({a,X q }).B q ], andP(a,X)Q 

where X = X' p n X' q and, for each i G {p, g}, X[ is the site Xj w.r.t. a global 
coordinate system, i.e. X[ = global(Xj,72-(5'i)). As stand-alone processes, 
P and Q can perform two compatible strong-split actions, namely p(a, X' p ) 
and p(a,X' q ) and evolve, resp., in S p [B p ] and S q [B q ]. As a consequence, 
P (a, X) Q becomes either S P [B P ] (a, X) Q or P (a, X) S p [B q ]. 

But these actions are compatible and, according to the intuition given 
so far, P and Q have to synchronize on their execution in order to split the 
bond (a, X). In other terms, a strong-split operation is enabled; such an op- 
eration must be performed before time can pass further and must produce as 
a result two independent 3D processes, i.e. the network of 3D processes (see 
Section [5J that contains both S p [B p ] and S q [B q ]. Similarly, we would allow 
synchronizations between compatible weak-split action in order to perform 
a weak-split operation. To properly deal with this kind of behaviours some 
technical details are still needed. We first allow synchronization on com- 

•ii t • i • i • • • i ■ p( a \X) u(a,X) 

patible split actions by introducing the transition relations =^ and => . 

Intuitively, we want that P(a,X)Q P ^ S P [B P ] (a, X) S q [B q ]. Now, we 
can 'physically' remove the bond (a, X) (this will be done by exploiting the 
function split over 3D processes we provide in the next section) and obtain 
the network of processes we are interested in. 

Definition 14 (Semantics of strong and weak splits) The SOS-rules 

that define the transition relations P =$■ 3DP x 3DP where p(a,X) E p(C) 
are given in Table ^ As usual, symmetric rules have been omitted. We 

also omit the rules defining the transition relations W 3DP x 3DP for 
u(a,X) G w(C) since these can be obtained from those in Table^by replacing 
each action p(-) with the corresponding action uj(-). 

Recall that strong-split operations require simultaneous split of multiple 
bonds. In this case, all the components involved in the reaction must all 
together be ready to synchronize on a proper set of compatible strong-split 
actions. Consider a more complex example 

P = S p [p({(a,X p ),(b,Y p )}).B p ], 

Q = S q [p({(a,X q )}).B q ], 

R = S r [p({(b,Y r )}).B r ], and (P (a, X) Q) (b,Y) R 
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where X = X' p n X' q and Y = Y p D Y^ (also in this case we write 
for i 6 {p, g, r}, to represent the site Xj in a global coordinate system). P 
can synchronize with Q and R to split, resp., the bonds (a,X) and (&,!''). 
Indeed, rules in Table [4] implies that 

(P(a,X) Q) (b,Y) R (S p [p({(b, Y p )}).B p ] (a, X) S q [B q ]) (b, Y) R 

P( ^ y) (S P [B P ] (a, X) S q [B q }) (b, Y) S r [B r ] 

After that, all 'pending strong-split requests' of (P (a, X) Q) (b, Y) R are 
satisfied. We say that such a compound process is able to complete a reaction. 
If it was, for instance, R = S r [e{t). P ({(b,Y r )}).B r ] then (P{a,X) Q) (b,Y) R 
would not have been able to complete a reaction, since (at least) one involved 
component, i.e. R, would not have been able to contribute to the reaction 
before than t units of time. In such a case, the bonds can not be spit at all. 
This concept is formalized by the definition below. 

Definition 15 (Bonds of 3DP-terms) The function bonds : 3DP — > p(C) 
returns the set of bonds that are currently established in P. It can be defined 
by induction on P £ 3DP: 

Basic: bonds(S[B\) = 

Comp: bonds(P (a, X) Q) = bonds(P) U bonds(Q) U {(a, X)} 

By an easy induction on P we can prove that P ^=^- ^ implies (a, X) € 
bonds(P). Moreover, we say that P E 3DP is able to complete a reaction, 

which we write as P \, iff either (1) P or (2) P P ^=> ^ Q for some 
p(a, X) and Q such that Q \ t . 



p p(^ pl q±^q. ae{a ,a } x = x p nx q 

StrSync — — 

P(a,X)Q p(a => } P' (a,X) Q' 

P P W ] P> 
StrPar — 

P(a,X)Q P % ' P' (a,X)Q 



Table 4: Transitional semantics for strong-split actions 
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Finally, if P is able to complete a reaction and there is at a least a bond 
that has to be strongly split (i.e. if P ^> does not hold), a reaction can 
actually take place and, as a consequence, time cannot pass further. Below 
we restrict the timed operational semantics of 3D processes as it has been 



defined in Def. 13 in order to take this aspect into account. 



Definition 16 Let P G 3D P. We say that P — > Q if P -w Q and either 
P ^> or P is not able to complete a reaction. 

Proposition [4] states that the function split is well-defined up to structural 
congruence over 3D processes we define below. 

Definition 17 (Structural congruence over 3D processes) We defi- 
ne the structural congruence over processes in 3DP 7 which we denote by =p, 
as the smallest relation that satisfies the following axioms: 

- S[B] = P S'[B] provided that S = s S' ; 

- P(a,X)Q= P Q(a,X)P; 

- P = P Q implies P (a, X) R = P Q {a, X) R; 

- Y C B(Q) n B(R) implies (P {a, X) Q) {b,Y)R= P P (a, X) (Q (6, Y) R). 

Proposition 4 Let P G 3DP well-formed. If (a,X) G bonds(P) there is a 
well-formed 3D process Q {a, X) R=p P. 

We also need the following closure result. 

Proposition 5 Let P, Q G 3DP and /j, G u(C) U p(C). Then: 

1. P => Q implies shape(Q) = shape(P); 

2. P well-formed and P Q implies Q well-formed. 

5 Networks or 3D processes 

Now we can define a network of 3D processes as a collection of 3D processes 
moving in the same 3D space. 

Definition 18 (Networks of 3D processes) The set N of networks of 
3D processes (3D networks, for short) is generated by the grammar: 

N ::= Nil I P I N \\N 
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where P G 3DP. Given a finite set of indexes I, we often write (|| Pi)ie/ 
to denote the network that consists of all Pi with i E I. We assume that 
I = implies (|| p) ie / = Nil. For N = ( || Pi) ie i we let Si = shape(P), 
for i G I, and define colliding(iV) as the set of all tuples (Pi,Pj,X) such 
that {Si,Sj,X) G colliding({5i} ie /) (see Defi A network N = (|| P) ie/ 
is said to be well-formed iff each Pi is well-formed and, for each pair of 
distinct processes Pi and Pj, the shapes Si and Sj do not interpenetrate. 
Moreover, we extend the definition of steer on networks in the natural way, 
i.e. steer t (|| Pi) ie i = (|| steer t P») i(E j. 

In our running example we construct a network of processes containing 
a proper number of HEX, ATP and GLC processes in order to replicate the 
conditions in a portion of cytoplasm. 

Definition 19 (Splitting bonds) The function split : 3DP x p(C) — > N 
is defined as follows: If {a,X} G bonds(P) n C and P =p Q {a,X} R then 
split(P, C) = split(Q,C) || split(P, C); if, otherwise, bonds(P) n C = 0, then 
split(P, C) = P. 

It is worth noting that split shapes maintain the same velocity until the 
next occurrence of a movement time step. As we mentioned above, this is 
not a problem because they will not trigger a collision and, thus, a shorter 
first time of contact. 

Proposition 6 Let P G 3DP well-formed and C C C. Then split(P, C) is a 

well-formed network of 3D processes. 

Definition 20 (Semantics of weak- and strong-split operation) Let 

P G 3DP a 3D process. If P \, we write that P A G N iff there is a non 
empty set of channels C = {(ai, Xi), ■ ■ ■ ,{a n ,X n )} C bonds(P) such that 

P = P p(a ^ Xl) P 1 ■ ■ ■ P{a ^ n) P n , P n andN = split(P n ,C7). Similarly, 

P -A iV iff there is a channel (a, X) G bonds(P) such that P =^ Q and 
N = sp\\t{Q,{(a,X)}). 

Since weak-split operations are due to a synchronization between just 
a pair of 3D processes, condition 'P is able to complete a reaction' is not 

needed, but 'P W ^ Q 1 suffices to our aim. 

Example 3 Let us consider P = H (atp, X) {A (glc, Y) G) where: 
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H = S h [p({(Btp,X ha ),(g\c,Y hg )}).HEX}, 

A = 5 [p({(atp,X oh )}).ADP + w(itp,X ah ).ATP], 

G = S g [p({(gk, Y gh )}).G6P + oj(gk, Y gh ).GLC], 

X' ha n X' , = X (here X' ha and X' , are the sites X^ a and X a h expressed 
w.r.t. a global coordinate system; this convention will be applied later on) 
and Y' hg n Y' gh = Y. According to the definitions given so far, P is able to 
complete a reaction since: 

P P( % X) S h [p({(g\c,Y hg )}).h\EX] (atp,X) (S^fADP] (glc,Y) G) 
p{g ^ Y) S h [HEX] (atp,X) (S a [ADP] (g\c,Y) S 9 [G6P}) = R 

Moreover, R ^> and 

split^C) = 5a[HEX]||split(5 a [ADP]<glc,y)5g[G6P],C) 
= S h [HEX]||(5 a [ADP]||5 fl [G6P]) = iV 

where C = {(atp,X), (glc, Y)}, implies P A N. Moreover, let us note that, 

for each t e T, P A P but P ^ since since P is able to complete a reaction 
and P does not hold. 

Below we define the temporal and functional behaviour of 3D networks. 
We assume that such networks perform basic actions that belong to set 
{oj, p, k}, where lj and p denote, resp., weak- and strong- split operations as 
a unique action (at the network level we only see whether shape bonds can 
be split or not) and k represents system evolutions due to collision response 



(see Section 5.1). We also let elements of the set {oo,p} U T to be ranged 
over by u, 1/ , ■ ■ ■ . 

Definition 21 (Temporal and Functional Behaviour of N-terms) 

Rules in Table\Q defines the transition relations — >-C N x N for t £ T and 
Ac N x N for v € {lj, p}. As usual, symmetric rules have been omitted. 

A timed trace from N is a finite sequence of steps of the form N = 
N Nx ■ ■ ■ N n = M. We also write that N A M if there is a timed 

n 

trace N = N N x ■ ■ ■ N n = M such that t = £{14 | ^ € T}. 

i=0 
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EMPTY/ PAFU 



N 4 N' M -4 M' 



Nil A Nil N\\mAn'\\M' 



N 4 JV' 
Par„ 



Table 5: Temporal and functional behaviour of 3D networks 



Proposition 7 Let t G T 7 P G 3DP 7 iV G N 7 raft P and AT well-formed. 

1. P N and P A iV implies N G N well-formed. 

2. N — > M implies M well-formed; 

3. N M implies M well-formed. 
5.1 Collision response 

In this section we describe the semantics of collisions response. As already 
said, the notion of compatibility between channels (and, hence, processes) 
has been introduced to distinguish between elastic and inelastic collision. In 
particular, collisions among compatible processes are always inelastic. So, 

if P - ' P ^> P' and Q - ' 9 \ Q' , with (a,X p ) and (a,X q ) compatible, and 
P and Q collide in the non-empty site X = X p f] X q we get a compound 
process (P' {a,X) Q')|vJ where the velocity v is provided by Equation (3) 
in Def. [9] Vice versa, a collision between two incompatible processes P and 
Q is treated as an elastic one. After such a collision, P and Q (actually the 
processes we obtain by updating their velocities according to Equations (1) 
and (2) in Def. [9J will proceed independently to each other. 

To resolve collisions, we introduce two different kinds of reduction re- 

laUons over 3D networks, namely and where P,Q are 

3D processes and X is a surface of contact (see Table [6]) . Intuitively, if 

N " % M [N K ') i M), then M is the 3D network we obtain once 
an elastic (inelastic, resp.) collision between P and Q in the surface of con- 
tact X has been resolved. These reduction relations also use the structural 
congruence over 3D networks. 
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Definition 22 (Structural congruence over 3D networks) 

The structural congruence over terms in N, that we denote with =, is the 
smallest relation that satisfies the following axioms: 

- Nil || N = N, N || M = M\\N and N || (M \\ R) = (N \\ M) \\ R; 

- P || N = Q || N provided that P = P Q. 

Rule ELAS in Table [6] simply changes velocities of two colliding but incom- 
patible processes guided by Equations (1) and (2) in Def. [9j while rule INEL 
joins two compatible processes P and Q to obtain a compound process whose 
velocity is given by Equation (3) in Def. [9] Note that we force P and Q to 
synchronize on the execution of two compatible actions (a, X p ) and (a, X q ) 
before joining them. In rules ELAS= and INELAS= we also consider struc- 



tural congruence over nets of processes. In Def. 23 we collect together all the 
reduction-steps needed to solve collisions listed in a given set of collisions 
colliding(iV); clearly N is a generic 3D network. 

Definition 23 (Resolving collisions) Let N G N and (P, Q, X) a tuple 
zncolliding(iV). N {P ' Q ' X) } M if either P ~ Q and N (P ' Q ' X) > . M or P ^ Q 
and N {P ' Q ' X \ e M. 

Moreover, we write that N -> M if either colliding(iV) = and N = M 
or colliding(iV) ^ and there is a finite sequence of reduction steps N = 
^ (PuQuX^ ^ (P k ,Q k ,X k)) ^ = M ^ thaf 

1. (Pi,Qi,Xi) e colliding(A r j_i) for each i e [l,k]; 

2. colliding(A7 fe ) = 0. 

Let also note that, at any given time t, colliding(A^) can be obtained 
from the set of all the pairs of processes in N that are touching at that time. 
This set and hence colliding(./V) is surely finite and changes only when we 
resolve some inelastic collision (this is because, after an inelastic collision one 
or more binding sites can possibly become internal points of a compound 
process, and hence are not available anymore). Moreover a collision between 
pairs of processes with the same shape can not be resolved twice. This 
is either because two processes P and Q have been bond in a compound 
process as a consequence of an inelastic collision, or because P and Q collide 
elastically and their velocities have been changed according to Equations 
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P A e Q = (Vp, V q ) 

ELAS (poY) 

(P || Q) || at S£^% e ( p|[ Vp j [|g[ v J)||JV 

N = N' N' (P,Q,X) ) e M 
ELAS= 



INEL- 



pi^p Q^^Q' a€{a,a} p X ^\Q = ^ 
(P || Q) || at < p ^ n ^ . ((P' (a, X p n X g ) Q')\[v]\) || iV 
JV = N' N' (p ' Q,x> > i M 

INEL= 



TV <^ M 



Table 6: Reaction rules for elastic and inelastic collisions 



(1) and (2) in Def. 
anymore (see Lemma 



in order to obtain two processes that do not collide 
in Appendix C). Thus, we can always decide if there 
is a finite sequence of reduction steps that allows us to resolve all collisions 
listed in colliding^) and hence obtain a network M with colliding(M) = 0. 



Proposition 8 Let N € N, P,Q G 3DP and X a not-emptyset subset o/P. 
Then N well-formed and N — — — > M implies M well-formed. 

By iterative applications of Proposition [8] (see Appendix C for the proof) 
it is also: 

Lemma 1 Let N a well-formed 3D network. Then TV A- M implies M 
well-formed. 

We are now ready to define how a network of 3D processes evolves by 
performing an infinite number of movement time steps. 



Definition 24 (System evolution) Let N,M G N and t,t' G T. We say 
that (N, t) =>• (M, t + r) iff one of the following conditions holds: 

1. t' = Ftoc(A0 < A and N 4 N' A N" and M = steer (t + t') N" ; 
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2. t' = A < Ftoc(iV) and N 4 N' and M = steer (t + if) N'. 
A system evolution is any infinite sequence of time steps of the form: 

i— 1 i 

(No, 0) 4 (N x ,ti) 4 (N 2 , h+t 2 ) — (Ni- 1} J2 tj) 4 (Ni, tj) ^ • • • 

j=l j=i 

Note that, for each i > 1, U = min{Ftoc(A^j_i), A} as discussed in 
Section [2j Moreover, in order to make sure that processes will never inter- 
penetrate during a system evolution, if Ftoc(iVj_i) < A, we first resolve all 
the collisions that happen after time U = Ftoc(A^-i) (by means of transition 
— >) and then apply the changes suggested by the function steer as described 
in Section [3.11 

Example 4 This example shows a possible evolution of the 3D network 
N = {HEX || ATP) || GLC where HEX, ATP and GLC are the 3D pro- 
cesses of Example [S| Below we use the following notation: 

- H(t) = (S h + t)[HEX\, A(t) = (S a + t)[NTP] and G(t) = (S g + t)[GLC] 
for each t G T. Note that HEX = H(0), ATP = A(0) and GLC = 
G(0); 

- C = p({(atp,X ha ),(g\c,X hg )}).HEX and, for any t < t h , HA(f) = 
u(atp, X ha ).HEX + e(t h - t).(g\c,X hg ).C; 

- AH(t) = w(atp,Xah).ATP + e(t a - t).p({Itp,X a/l }).ADP for anyteT 
with t < t a ; 

- GH(t) = w(gk,X 5h ).GLC + e(t g - t).p({gk, X gh }).G<5P for any t < t g . 

Let t\ = Ftoc(A?o) and assume t\ < A. By the operational rules, it is 
N 4 H(tx)\\A(ti)\\G(ti) = Nt). We also assume that colliding(iVQ) = 
{(fl-(ti),A(ti),Jf)} where X = X' ha C\X' ah / 0, X' ha = g\obz\(X ha ,K{S h + 
h)) and X' ah = g\obz\{X ah ,K{S a + t x )). Then: N' Q 4 P(h) \\ G(h) = N(f 
where P(h) = ((S h + ti)[HA(0)](atp, X)(S a + ti)[AH(0)]) [ Vfeo ]| and w ha = 

H(t\) 4^ A(ti). Finally: 

(N ,0) ^ (N u ti) 

where N x = steer t\ Nq = steer t\ P(t\) || steer t\ G(t\) = 
^(*i)|[vi]| ||G(ti)|[v 2 ]|. Note that: 
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P(ti)\[vi}\ = ((5 ft + ti)[HA(0)](atp,X)(5 a + t 1 )[AH(0)])[vil 
= (^[HA(0)](atp,X)^[AH(0)]) 

where Si = ((S A +ti)[vil) and = ((S +ti)|[vi]). Moreover, G(ti)|[v 2 ] = 
((5 s + i 1 )|[v 2 |[GLC] = ^ 1 [GLC]. 

Let t2 = Ftoc(iVi) and assume t 2 = < min{t a , Ajp^j Below we write 
G'(t 2 ) and P'(t2) to denote, respectively, the 3D processes (S^+t^fGLC] and 
[(Sl+t 2 )[HA(t h )](atp,X+t 2 -v 1 )(Sl+t 2 )[AH(t 2 )]). Again by the operational 

rules, Ni % P'(t 2 ) \\ G'(t 2 ) = N[. 

Let colliding^) = {(P'(t 2 ),G'(t 2 ),Y)} where Y = X' hg n X' gh + 0, 
X' hg = g\oba\(X hg ,n(Sj l + t 2 )) andX gh = g\ote\{X gh ,K{Sl + t 2 )) C 

B(P'(t 2 )). Ifv gh = P'(t 2 ) A< G'(t 2 ), then 

N[ A (P(t 2 ) (glc,K) G(t 2 ))[v flfc ] = iV{' w/iere G(t 2 ) = (S l g + t 2 )[GH(0)] and 
P(t 2 ) = ((5^ + t2)[C]<atp,X + t 2 -vi)(5i+t 2 )[AH(t 2 )]). fina% 

(JVi,ti) 4 (N 2 ,tt + t 2 ) 

where N 2 = steer (t x + i 2 ) = [P(t 2 ) (glc,F) G(t 2 )) |[v 3 ]|. Observe that: 
(P(t 2 ) (glc,F)G(i2))|[v 3 ]| = (P(t 2 )|[v 3 |) (glc,y) (G(t 2 )|[v 3 ]|) = 
((S 2 [C] (atp, X + t 2 • Vl ) S 2 [AH(i 2 )]) (glc, y) S 2 [GH(0)] 

i^ere = (S* + i 2 )|[v 3 ], S 2 = + t 2 )[v 3 | and S 2 g = + t 2 )|[v 3 |. 

At this stage the network contains just one process and, as a consequence, 

no collisions are possible. Thus, Ftoc(iV 2 ) = oo. Assume t g = t a — 1 2 < A 14 

If we let X g = {X + t 2 ■ vi) + t g ■ v 3 and Y g = Y + t g ■ v 3 , then: 

N 2 % ((S\ + t g )[C] (atp,X g ) (S 2 a + t g )[AH(t a )}) (g\c,Y g )(S 2 g + t g )[GH(t g )] A 

((S 2 h + t g )[HEX] || (S 2 + t 9 )[ADP]) || (Si + t,)[G6P] ^ 

((Si + A) [HEX] || (S 2 + A)[ADP]) || (S 2 + A)[G6P] = N 2 . Thus: 

(N 2 ,t 1 + t 2 )^(N 3 ,t 1 +t2 + A) 

13 If were ti < th- the 3D processes HA(t2) and GLC would be no more compatible, and 
a collision between them would be treated as elastic. On the other hand, if were £2 = t a 
the idling time for AH(t a ) would be over. As a consequence, time would pass further only 
after the execution of a weak-split operation that splits the bond between the Hexokinases 
and the Atp molecules 

14 If were t g 7^ t a —ti the reaction could never proceed since the involved molecules would 
never be able to release - all together - the bonds. Thus the system would deadlock. 
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where N 3 = steer fa + t 2 + A) N' 2 = (Sjj[HEX] || ^[ADP]) || 5|[G6P], 
Sf = (5f+A)|[vi| and V; = steer fa+i 2 + A) (Sf + A) for each i G {h,a,g}. 



We can prove the following basic property of the Shape Calculus stat- 
ing that any system evolution does not introduce space inconsistencies like 
interpenetration of 3D processes or not well-formed processes. 



Theorem 1 (Closure w.r.t. well-formedness) Let N be a well-formed 

ft 

network of 3D processes. If (N, t) =4> (M, t + t') then M is well-formed. 
Proof: Assume that (N, t) 4 (M, t + t') because of JV 4 N' A N" and 



M = steer (t + t') N" (the other case - see Def. 24 - is similar). Then, by 



Proposition [7] (item 3) and Lemma [lj well-formed and A^ 
implies A^" and hence M = steer (t + 1') N" well-formed. 



4 N' A N" 



6 Conclusions and Future Work 

We have defined the full timed operational semantics of the Shape Calculus 
and we have introduced a notion of well-formedness of the different objects 
of the calculus. We proved that the evolution of a well-formed network of 
3D processes is always well-formed, that is to say, no spatial or temporal in- 
consistencies can be introduced by the dynamics of the calculus. As future 
work we intend to provide verification techniques for the Shape Calculus. 
In order to do this we believe that a sort of "tailoring" should be made on 
the calculus, making some parts (e.g. movement) more abstract and other 
parts (e.g. behaviours) more specific adding quantitative information (for 
instance probabilities or costs). The whole process will then be supported 
by the definition of proper logical languages to specify properties of inter- 
ests. Of course we expect that some approximations will be necessary to 
reach computability and/or feasibility. Another direction of future work 
is the possibility to include in the calculus some new useful, but in some 
cases complex, concepts such as re-shaping, message passing of values, and 
communication by perception of a compatible process in the neighbourhood. 
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Appendix A: Proofs of Section [3] 

The following lemmas are first needed. 

Lemma 2 Let S G § well-formed and t G T. Then: 
1. V (S + t) = v(5) and m(5 + 1) = m(,S); 
S. 71(5 + 1) = Tl(S) + t • v(S); 

3. V(S + t) = V{S) + t • v(5) fand, Zience, £(S + t) = B(S) + t • v(S) ). 

Proof: We prove Items 1, 2 and 3 by induction on S G § well-formed. 

Basic: S = (V, m, p, v), S + t = (V + (t ■ v), m, p + (t • v), v) and: 

1. v(5 + i) = v = v(S) and m(S + t) = m = m(S); 

2. K(S + t) = p + (t • v) = ft(S) + (t • v); 

3. 7?(S + t) = F + (t • v(5)) = + (t • v(5)). 

Comp: S = St (X) S 2 and 5 + 1 = (Si + 1) (Y) (S 2 + t) where Y = X + 1 ■ v 
and v = v(5) = v(Si). By induction hypothesis: 

1. v(Si + i) = v(5i) = v and m(5j + t) = m(Si) for i = 1,2. Thus: 
v(5-|-i) = v and m(S+t) = m(Si+i)+m(S2+t) = m(Si)+m(S2) = m(S). 

E m(Si+t)-Tl(Si+t) E m(Si)-(n(Si)+t-\r) T, m(S l )-n(S l ) 

2. n{s + t) = i -^— 2 = ^ — - 2 = ^4 + 

E m{Si+i) E m(Si) E m(Si) 

i. — 1 i — 1 z — 1 

t • v = K(S) + t • v 

3. V(S + t) = P(5i + 1) U P(S 2 + 1) = 7>(5i) U 7>(S 2 ) = 

To prove Prop. [T] we also need the following lemma; its proof has been 
omitted because it is similar to that of Lemma [2] 

Lemma 3 Let S G § well-formed and w € V. Then: 

1. v(5|[w]|) = {w} and m(5|[w]) = m(5); 

2. P(5'|[w]|) = 7>(S) (and, hence, B(5|[w]) = B(S)). 

Proposition [T] Lei 5 G S, f G T and w G V. If S is well-formed then S + t 
and S\[wJ are well-formed. 

Proof: We prove these statement by induction in S G §. 

Basic: S = (V, m, p,v). Both S + t = (V + t ■ v, m, p + t ■ v, v) and 
5|[w| = (V, m, p, w) are well-formed shapes. 

Comp: S = St{X)S 2 . Let v = v(S) = v(5<) for i = 1,2. By Defs J5| 
and§ it is 5 + t = (Si + t) (Y) (S 2 + t) where Y = X + t ■ v, and S\[w]\ = 
(<Si|[w|) (X) (^([w]!). We first prove S + t wellformedness. 
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1. By induction hypothesis, S\ + t and S 2 + 1 are well-formed; 

2. S well-formed implies that X is non-empty set that is equal to V{S\) fl 
V{S 2 ) = B(Si) n B(S 2 ). By Lemma [2j we also have that Y = X + t • v is 
non-empty and F = V(Si + 1) D P(S 2 + t) = B(5i + t) D £(S 2 + i). 

3. By Lemma [2j-l, v(iS + t) = v where v = v(S*j + t) = v(S'j) for i = 1,2. 

Similarly we can prove that S\[w]\ is well-formed. Shapes 5i|[w| and S^w] 
are well-formed (this follows by induction hypothesis); by Lemma [3]- 1 it is 
v(S\[w}\) = w; finally ^ X = V{Si) n T(S 2 ) = B(Si) D B(S 2 ), V(Si) = 
V(Si\[w]\) and B(Si) = B(Si\[w]\) for i = 1,2 (by Lenuna§2) imply / 

x = P(Si\[w]) n v(S 2 \[w]\) = B(5i|[w]|) n B(5 2 [wJ). 
Appendix B: Proofs of Section [4] 

This section is devoted to prove main results stated in Section |4j 

Proposition [3] Let P £ 3DP well-formed. Either P Q or P Q implies 
Q £ 3DP well-formed. 

Proof: We first prove (by induction on P E 3DP) that: 

1. P Q implies shape(Q) = shape(P) + 1; 

2. P AQ implies shape(Q) = shape(P). 

Basic: P = S[B\. 

1. By operational rules, P A Q implies P> A P/ and Q = (S + i) [-B'] . Thus: 
shape(Q) = S + t = shape(P) + t. 

2. P A Q implies Q = S[B'} for a proper 6' £ B (see rule BASIC C in 
Table g). Hence shape(Q) = 5 = shape(P). 

Comp: P = Pi {a,X)P 2 . 

1. P A Q implies Pi A Qi, P 2 A Q 2 , and Q = Qi (a, F) Q 2 where y = 

(i • v(P)). By induction hypothesis, shape(Q) = shape(Qi) (Y) shape(Q 2 ) = 
(shape(Pi) + t) (Y) (shape(P 2 ) +t) = shape(P) + t. 

2. P A Q implies either P 1 A Q x and Q = Q 1 (a, X) P 2 or P 2 A 
Q 2 and Q = P\(a,X)Q 2 . Let us consider the former case (the latter 
one can be proved similarly). By induction hypothesis it is shape(Q) = 
shape(Qi) (X) shape(P 2 ) = shape(Pi) (X) shape(P 2 ) = shape(P). 

Now we prove that P well-formed and P A Q imply Q well-formed. Again, 
we proceed by induction on P well-formed. 
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Basic: P = S[B}. If P ^ Q then Q = (S + t)[B'\ where P A In 
this case, S well-formed and Prop. [T] imply shape(Q) = S + 1 well-formed. 
Moreover global(X, 1Z(P)) C B(P) for each site ICP that occurs in B and, 
hence, in B' . By Lemma|2}4, global(X, Tl{Q)) = global(X,ft(P)+i-v(P)) = 
global(X,ft(P))+i-v(P) C B{P)+t-v(P) = B(Q). Thus, Q is well-formed. 

Comp: P = Pi (a, X) P 2 with P well-formed for i = 1, 2 and X C B(Pi) n 
£(P 2 ). P A Q implies P A Q i5 for j = 1,2, and Q = Qi (a,Y) Q 2 with 
Y = X + (t ■ v(P)). By induction hypothesis, Q\ and Q 2 are well-formed. 
Moreover, X C £(Pi)n£(P 2 ) implies Y = X + t- v(P) C (fi(Pi) n B{P 2 )) + 
t ■ v(P) = (B(Py) + t • v(P)) n (£(P 2 ) + t • v(P)) = B(Qi) n B{Q 2 ). 

It remains to prove, by induction on P, that P well-formed and P A Q 
imply Q well-formed. 

Paszc: P = S[B\. If P A Q then P A P' for a proper P' € 3DP and 
Q = S[B']. P well-formed implies S well-formed. Moreover, for each X that 
occurs in B' (and hence in P), global(X, K(P)) C £(P) = B(£) = B(Q). 
So, Q is well-formed. 

Comp: P = Pi (a, X) P 2 with P well-formed for % = 1, 2 and X C B(Pi) n 
S(P 2 ). P A Q implies either Pi A Qi and Q = Qi (a,X) P 2 or P 2 A Q 2 
and Q = P\ (a, X) Q 2 . We only prove the former case (the latter one is 
similar). By induction hypothesis Q\ is well-formed; moreover, shape(Qi) = 
shape(Pi) and, hence, B{Qi) = B{Pi). Thus, Q\ and P 2 are well-formed 
with X C B(Px) n B(P 2 ) = B(Qi) n B(P 2 ), i.e. Q is well-formed. 

Proposition [4] Let P £ 3DP well-formed. If (a,X) G bonds(P) i/iere is a 
well-formed 3D process Q (a, X) R=p P. 

Proof: By induction on P well-formed. 

Basic: P = S[B]. This case is not possible since bonds(P) = 0. 

Comp: P = Pi (6, Y) P 2 with P well-formed for i = 1, 2, Y C 0(Pi) n£(P 2 ) 
and v(Pi) = v(P 2 ) = v(P) = v. If (a,X) = (b,Y) it suffices to choose 
Q = Pi and P = P 2 . Assume (a, X) G bonds(Pi) (the case in which (a, X) G 
bonds(P 2 ) is symmetric). By induction hypothesis, Q\ (a,X) R± =p P\ is 
well-formed and P = Pi (b, Y)P 2 = P (Q x (a, X) Pi) (b, Y) P 2 . We have two 
possible subcases: 

1. YQB(Ri) and P = P Q (a,X) P where Q = Q x and R = R\ (b, Y) P 2 , 

2. Y C B(Qi) and P = P Q {a, X) R where Q = R ± and P = Qi {b, Y) P 2 . 
We only consider the former case (the latter one is similar) and prove that 
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Q{a,X)R is well-formed. To this aim we have to show that: (i) R = 
R\ (b, Y) P2 is well-formed (note that Pi =p Q\ {a,X)R\ well-formed and 
Q = Qi implies Q\ well-formed), (it) v(Q) = v(P) (this follows easily 
because P well-formed implies v(Qi) (i.e. v(Q))= v(Pi) = v(P 2 ) = V (R)) 
and (Hi) X C B(Q) n B(R). Below we prove first (z) and then (ra) 

Pl =p Ql (a>, X) R\ and P = P± {b, Y) P2 well-formed imply, resp., R± and 
Pi well-formed; moreover, v(R\) = v(Pi) = v = v(P 2 ) (because P\ =p 
Qi (a, X) R\ is well- formed) . To show that R = R\ {b, Y) P2 is well-formed, 
it remains to prove that Y C B(R\) D 0(P 2 ). We proceed by contradiction. 
Assume y G Y such that y £ £(i?i) n B(P 2 ). Then, F C B(R 1 ) implies 
y ^ B(P2) and, hence, y ^ B(P\) n ,6(P 2 ). But, this is impossible because 

rc n£(P 2 ). 

We prove that X C £(Q) n B(R) = B(Q{) n B(P) by contradiction. Recall 
that Qi (a,X) R\ well-formed implies X C B(Qi) C\B(R\). Assume that 
there is x G X C B(Qi) n B(Pi) such that x ^ n B(R). Then 

x £ B(R) = {B{R 1 )UB(P 2 ))\{z\z is interior of V(Ri) UP(P 2 )}. Moreover, 
x G X C S(Pi) C B(P X ) n 0(P 2 ) and x £ 0(P) implies that x is an 
interior point of V(R\) U P(P 2 ), i.e. x is an interior point of P(P 2 ); this 
is because x G B(R\) implies that x cannot be an interior point of V(R\). 
Finally, x G X implies that x is also an interior point of Pi =p Q\ {a, X) R\. 
Thus, Pi and P 2 interpenetrate each other. But this is not possible since 
P = P Pi (b, Y) P 2 is well-formed. 

Proposition ^pet P,Q G 3DP and p G uj(C) U p(C). Then: 

1. P Q implies shape(Q) = shape(P); 

2. P well-formed and P Q implies Q well- formed. 

Proof: We only prove the statement for p = p(a,X); if p = uj(a,X) the 
statement can be proved similarly. 

Basic: P = S[B]. This case is not possible since P f^%- \ 

Comp: P = Pi (b, Y) P 2 with P well-formed for i = 1, 2 and X C B(Pi) D 
£>(P 2 ). We distinguish two possible subcases: 

{1) (b,Y) = (a,X), Pi p{a ' Xl \ Q 1; p 2 p(g ' Xa) > Q 2 with a G {a, a} and 
X = X x n X 2 , and = Qi (a, X) Q 2 . 

1. By Prop. [3^-2 , shape(Qi) = shape(Pj) for i = 1,2 and, hence, shape(Q) = 
shape(Qi) (X) shape(Q 2 ) = shape(Pi) (X) shape(P 2 ) = shape(P). 
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2. By induction hypothesis, Pj well-formed implies Qi well-formed, for i = 
1,2. Moreover X C 0(P X ) n 0(P 2 ) and, again, shape(Qj) = shape(P) for 
i = 1, 2 imply X C £?(Qi) n B(Q 2 ). So the 3D process Q is well-formed. 

(m) Either Pi p( ^° Qi and Q = Qi (a, X) P 2 or P 2 p( ^° Q 2 and Q = 
Pi(a,X) Q 2 . 

1. By induction hypothesis it is sbape(Qi) = shape(Pj) for i = 1,2. As in 
the previous case, we can prove that shape(Q) = shape(P). 

2. By induction hypothesis we have that Qi and Pj, for {i,j} £ {1)2}, 
are well-formed. Moreover, by Item 1, shape(Qj) = shape(Pj), and hence 
X C B(P 1 )nB(P 2 ) = B(Qi)C]B(Pj). Also in this case we can conclude that 
Q is well-formed. 

Appendix C: Proofs of Section [5] 

Proposition [6] Let P € 3DP well-formed and C QC. Then split(P, C) is a 

well-formed network of 3D processes . 

Proof: By induction on the number of channels in bonds(P) Pi C. 

bonds(P) n C = 0. In such a case split(P, C) = P is well-formed. 

(a, X) e bonds(P) n C / 0. By Prop. [3] there are Pi , P 2 well-formed 
such that P = P Pi (a,X)P 2 and split(P,C) = split(Pi,C) || split(P 2 , C) = 
N\ || N 2 . By induction hypothesis iVj is a well-formed 3D network for i = 1,2. 
Now, in order to prove that split(P, C) is well-formed, it still remains to 
prove that if Q\ and Q 2 are two 3D processes in the network N\ and N 2 , 
respectively, then Q\ and Q 2 do not interpenetrate. Assume, towards a 
contradiction, that there are Q\ and Q 2 in N% and N 2 , respectively, such 
that Qi and Q 2 interpenetrate, i.e. such that there is (at least) a point x 
that is interior of both Q\ and Q 2 . 

For each i = 1,2, it holds that either bonds(Pj) n C = 0, N { = Pi 
and Qi = Pi or bonds(Pj) n C = { {a x ,Xi) , {a 2 ,X 2 } , .., (a ni ,X ni ) }, P { = P 
P/ {at, Xi) (Pf (a 2 , X 2 ) ..(P™ (a ni , X ra ) P™ +1 )) (this follows by iterating 
Prop. Q, iVj = P/ || (Pf || • • • (P™ || P" i+1 )) and = P/ for some j = 
1, 2, • • • , ni,ni + 1. Thus, in both cases the set of interior points of Qi is a 
subset of interior point of Pj. As a consequence, if Q\ and Q 2 interpenetrate 
so do Pi and P 2 . But this is not possible because P is well-formed. To 
prove Prop. [7] we need the following preliminary result. 
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Lemma 4 Let N,M G N, P,Q G 3DP arwi t e T suc/i i/iai TV -4 M and 
P 4 Q. The network N contains P iff M contains Q. 

Proof: By induction on N G N. 

iV = Nil. iV 4 Nil and Nil does not contain any 3D process. 

N = P with P G 3DP. By our operational semantics, N 4 M iff M = Q 
and P -> Q. 

iV = JVt || jV 2 . iV 4 Mi || M 2 = M iff JV 4 4 Mi for i = 1,2. P is contained 
in iV iff P is contained either in N± or in N 2 iff (by induction hypothesis) 
and Q is contained either in M\ or in M2, i.e. iff Q is contained in M. 

Proposition [7] Let £ G T, P G 3DP, AT G N, with P and N well-formed. 

1. P N and P A N implies N G N well-formed. 

2. N 4 M implies M well-formed; 
t 



3. N =>■ M implies M well-formed. 



Proof: Item 1 can be proved by (iterative applications of) Prop. [5]-2 and 



Prop. [6] Item 3 follows directly from Items 1 and 2 (see Definition 21). In 
what follows, we prove Item 2 by induction on N G N. 

N = Nil. In such a case iV 4 M implies M = Nil that is well-formed. 

N = P with P G 3DP. The statement is a trivial consequence of Prop. |3]-3. 

JV = JVi || N 2 . If JV 4 M then 4 Mj, for t = 1, 2, and M = Mi || M 2 . By 
induction hypothesis, Mi and M 2 are well-formed. So, it remains to prove 
that if Qi and Q2 are 3D processes that compose the networks Mi and M 2 , 
resp., then Q\ and Q 2 do not interpenetrate. We proceed by contradiction. 
Assume that there are Q± in Mi and Q2 in M 2 that interpenetrate each 
other. By Lemma |4j Mj contains Qi iff iVj_ contains P, for some P, G 3DP 
such that Pj — > Qj. Moreover, by Prop, pl-l, shape(<3i) = shape(Pj) + t. 
Thus, if Qi and Q 2 interpenetrate the same do Pi and P 2 . But this not 
possible since ./V is well-formed. 

Lemma 5 Let N, M G N such that N is well-formed and N — — — > M for 
some (P,Q,X) G colliding(iV). If P' and Q' are two processes composing M 
with shape(P') = shape(P) and shape(Q') = shape(<3), then (P',Q',Y) £ 
colliding(M) for any 7 CP. 
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(P Q X) (p Q x) 

Proof: By rules in Table 6 N — — — > e M implies N — — — > e M implies 

M ee (P|[v p ]| || Q\[v q ]\) || N' where the pair of velocity (v p ,v g ) = P Q- 
Moreover, Equations (1) and (2) in Def.[9]ensure that -P|[v p | and Q|[v g | (since 

N is well-formed, these are the only processes in M with the same shape of 

(_P Q X) 

P and Q, resp.) can not collide anymore. Now assume TV" — — — >i M. By 



Def. 



23 



there are compatible (a, X p ) , (a, X q ) G C such that P - ' P \ P', 



Q Q' and M = {{P' (a,X) Q')\[v]\) \\ N' with a G {a, a}, X = X p n 

X q ^ and v = P i — >{ Q. The statement follows because there are no P 
and Q' in M such that shape(P') = shape(P) and shape(Q') = shape(Q). 

Proposition [8] Let JVgN, P, Q G 3DP anc? X a not-emptyset subset off. 

(P Q X) 

Then N well-formed and N — — — > M implies M well-formed. 

Proof: If N {P ' Q ' X \ e m then N = (P \\ Q) \\ N' , M = (P\[v p }\ \\ Q\[v q }\) \\ N' 

and (v p , v g ) = P < — > e Q. Moreover, N well-formed implies that P, Q and 
N' are well-formed and processes P and Q do not interpenetrate each other 



and with any other 3D process composing N' (see Def. 18). Finally: -P|[v p ] 
and Q|[v 9 ]| are well-formed 3D processes with shape(i- , |[v p |) = shape(P) and 
shape(Q|[vg|) = shape(Q) (see Lemma [3]). Thus, P|[v p | and Q|[v 9 ] do not 
interpenetrate each other and with any other 3D process in N' and we can 
conclude that M is well-formed. 

(P Q x) 

Now assume gain, N = (P \\ Q) \\ N' and, by rule INEL, 

there are (a, X p ), (a, X q ) G C compatible such that a G {a, a}, P - ' F \ P', 
Q Q' and M = ((P 1 (a,X) Q')[vJ) || N' where X = X p D X q ^ $ 

and v = P < X > i Q. By Prop. [s]-3, P and Q well-formed implies P' and Q' 
well-formed. Moreover: 

(1) P P' and Q Q' imply X p C B(P) and X q C B(Q), resp. 
Thus, I = I„nl, C B{P) n B(Q) = B(P') n B(Q') (since Prop. [3}2 
implies shape(P') = shape(P) and shape(Q') = shape(Q)), P' {a,X) Q' and 
{P' {a,X) Q')l[ v ]l are well-formed 3D processes. 

(2) shape((P' (a, X) Q')\[v}\) = shape(P' (a, X) Q') = V(P') U V(Q') = 
V{P) UV(Q). As a consequence, (P' (a, X) Q')\[v]\ can not interpenetrate 
any 3D process that compose the network N' (otherwise either P or Q must 
do the same). 

Aso in this case we can conclude that M is well-formed. 



40 



